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o3 ' Abstract 

We study bubbling geometry in topological string theory. Specifically, we analyse Chern- 
Simons theory on both the 3-sphere and lens spaces in the presence of a Wilson loop of 
an arbitrary representation. For each three manifold, we formulate a multi-matrix model 
whose partition function is the Wilson loop vev and compute the spectral curve. This 
spectral curve is closely related to the Calabi-Yau threefold which is the gravitational 
dual of the Wilson loop. Namely, it is the reduction to two dimensions of the mirror to 
the Calabi-Yau. For lens spaces the dual geometries are new. We comment on a similar 
matrix model relevant for Wilson loops in AdS/CFT. 
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1 Introduction and summary 

A useful aspect of duality between a gauge theory and a gravitational system is the 
emergence of spacetime through dynamics of gauge theory. Deeper understanding of 
emergent geometry should help us find new formulations of string theory and quantum 
gravity that may be used to address fundamental questions in physics. 

In gauge/gravity duality, the vacuum state corresponds to a certain background space- 
time, and inserted operators to excitations. The fields of gauge theory backreact sig- 
nificantly to the insertion of some operators. The corresponding gravitational dual is a 

1 



new geometry that shares the asymptotics with the original background. A bubble of 
new cycles supported by flux appears, and the new spacetime is thus called the bubbling 
geometry. The bubbling phenomenon was originally found for local operators [1], and 
was generalized to Wilson loops [21 El H] in AdS/CFT. It is useful to introduce a ma- 
trix model which captures the dynamics of all the relevant fields that respond to the 
operator insertion [71 [8] . One is able to visualize the backreaction in terms of eigenvalue 
distributions, which in turn encode the bubbling geometry on the gravity side. 

The current work studies the topological string version of bubbling phenomena [5], 
which naturally extend the Gopakumar-Vafa gauge/gravity duality [6] . More specifically 
we consider U{N) Chern-Simons theory on S^ or lens space L{p, 1) = S^/Zp with Wilson 
loop insertions. The Wilson loop operator is defined as 



W, 



R 
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where A is the gauge field and is integrated along the unknot. For S^fZp we take the 
unknot that generates the fundamental group. The trace is evaluated in an arbitrary 
representation R of U{N). Throughout the paper the symbol R also denotes the corre- 
sponding Young tableau, and we parametrize it as in Figure [H Each edge length be it 
nj or kj, will correspond to the size of a new cycle in the bubbling geometry. 




Figure 1: The Young tableau R, shown rotated and inverted, is specified by the lengths 
rij and kj of the edges. Equivalently, rij and kj denote the lengths of the black and 
white regions that are obtained by vertically projecting down the edges in R onto the 
horizontal line, rim+i is defined by YlT=i "^i ~ ^ ■ 



Building on the earlier work [9|, [TO] , we formulate a matrix model whose partition func- 
tion is the vev of the Wilson loop in S'^ or S'^ /Xp. We then study the eigenvalue dynamics 
in the large A^ limit and derive the spectral curve. For S^ the spectral curve is precisely 
the mirror of the bubbling toric Calabi-Yau geometry identified as the gravitational dual 
of the Wilson loop in [5] . The topology of this threefold depends on the data encoded in 
the Young tableau R: its toric web diagram is shown in Figure W^a) 



For the lens spaces S'^/Zp, the backreaction of the fields to the Wilson loop leads to 
additional classical vacua, and the path-integral splits into sectors corresponding to the 
different vacua. Because the matrix model we formulate computes the Wilson loop vev 
in each sector, we propose that for given N,p, and i?, a single Wilson loop insertion is 
dual to a sum over bubbling geometries. Each term in the sum is the toric Calabi-Yau 
that is mirror to the spectral curve which we derive. The summed geometries have the 
same toric data shown in Figure [2](b)lj except different values of Kahler moduli. As in 
the S'^ case, the topology of the geometry depends on the Young tableau data. 






(a) (b) 

Figure 2: (a) The toric web diagram for the bubbling Calabi-Yau dual to the Wilson loop 
Wr in 5*^. It has 2m + 1 copies of P^. (b) The web diagram for the bubbling Calabi-Yau 
dual to Wr in lens space S'^/Zp with p = 3. The diagram is a chain of m + 1 basic units. 

The paper is organized as follows. Section [2] focuses on the S*^ case. In subsection 12.11 
we present the matrix model for a Wilson loop in S^. In subsection 12.21 we derive the 
matrix model from physical arguments. Specifically we present it as an open string field 
theory of a D-brane configuration that realizes the Wilson loop. Then we algebraically 
derive the matrix model in subsection 12.31 In subsection 12.41 we solve the matrix model 
in the large N limit and derive the spectral curve, which is the mirror of the bubbling 
Calabi-Yau found in [5]. 

Section [3] deals with lens space S'^/Zp, and is structured in parallel with section [2l 
For each vacuum of the gauge theory with Wilson loop insertion, we derive the spectral 
curve. We propose that the mirror toric Calabi-Yau is the bubbling geometry dual to 
the Wilson loop. 

Appendix |A] summarizes the notation regarding the Young tableau data. In appendix 
Owe study alternative matrix models that compute the Wilson loop vev. The models are 
the direct analog of the matrix models for A/" = 4 Yang-Mills considered in [11] . Appendix 
[P] is targeted at readers interested in AdS/CFT. We use the algebraic techniques in 
subsection 12.31 to formulate a matrix model, whose partition function is the vev of the 
supersymmetric circular Wilson loop in A/" = 4 Yang-Mills. In this formulation it is very 
easy to derive the eigenvalue distributions for the Wilson loop found in [121 [131 El E]- 



^ To take the limit p ^ 1 in Figure [2l one need to apply an SL{2,'Z) transformation. 



2 Bubbling Calabi-Yau for S^ from a matrix model 

2.1 Matrix model for a Wilson loop in S^ 

The realization that the open topological A-model can be reduced to a matrix model 
first appeared in Marino's work [9], and a B-model version of this idea was subsequently 
derived by Dijkgraaf and Vafa [ISJ. Both derivations are of course mirror to each other 
as was demonstrated for certain examples in the nice work [10]. We are interested here 
in the A-model, which is of course equivalent to Chern-Simons theory |16], possibly with 
instanton corrections p^fTTj fT8]. 

Marino's observation for Chern-Simons theory on S^ with the gauge group G was that 
the partition function is 

Z = dnue ^s^ 



\,\{du..\[[2.u.h'^) e-^^'-s (2.2) 

i— 1 i^o V / 



where the topological string coupling constant Qs is identified with the Chern-Simons 
coupling constant, U = e'^, and dnu is the Haar measure on G with unusual integration 
range. On the second line we specialized to the case G = U{N), and each Ui is integrated 
from — oo to +oo. This observation by itself may not be overwhelming since it is a 
reformulation of Witten's classic result for the partition function [19]. The main utility is 
the generalization to different manifolds, where they carry topological data [9], [ini EO, EI] 
and to Wilson loops which we describe in the current work. 

One very interesting feature of fl2.2p however is that it secretly knows about the geo- 
metric transition of Gopakumar and Vafa [6] . While Chern-Simons theory is equivalent to 
the open A-model on the deformed conifold, the spectral curve of (12.21) is directly related 
to the resolved conifold. If the Calabi-Yau threefold mirror to the resolved conifold is 
defined by the equation xy = /(e", e^), the spectral curve is then given by /(e", e"") = 0. 
The orientifold case was worked out in [22]. Our main interest in this paper is to gener- 
alize this aspect to include the insertion of Wilson loop operators. Wilson loops in the 
topological gauge/gravity duality have been considered before by Ooguri and Vafa [23] . 
The current work and the previous work [5], [2ll [25] extends this in two ways. Firstly, 
the full backreaction of the Wilson loop is taken into account, as explained in [23] this 
means the Wilson loop vev can be expressed in terms of purely closed string enumerative 
invariants. Secondly we provide a dictionary for a single Wilson loop in a particular 
representation R, whereas in [23] a sum of Wilson loop insertions was considered where 
the summation is over representations. 



Wilson loop operators 

Wr = TrPe^^ (2.3) 

are specified by two pieces of data: the representation R of the gauge group G and a curve 
7 in M which the gauge field is integrated over. We will be considering all representations 
of U{N) such that the nj and kj in Figured] are large and our 7 will be the unknot. 

The relation between Chern-Simons theory and the matrix model was extended to 
include Wilson loops in p!0] : 

{Wr) = J dHue-^^'^'^'hiRe^ (2.4) 

We will show that the vev of the Wilson loop is in fact the partition function of the 
following matrix model: 

/m+l 
1=1 I<J 

(I) (/) \ 2 



nu^n^-s^n 2-h. 



^ ^-557E>(<')%^^E.«f 



/ ^(^)_^(-^)\ 

xJin 2sinh^^-^ , (2.5) 

KJ ij \ / 

with 

m -. 7—1 , m+l 

Li = ^kj--^nj + -^ nj for/ = l,...,m + l. (2.6) 

j=i J=i J=i+i 

This is a Gaussian {m + l)-matrix model with certain interactions which in the next 
section we explain from the target space viewpoint. 

2.2 Physical derivation of the matrix model 

In this subsection we derive the matrix model (12.51) as the world-volume theory in a 
D-brane configuration that is equivalent to the Wilson loop insertion. Further geometric 
transition of the branes leads to the purely closed string geometry in Figure [2t^a), and 
the three steps are summarized in Figure [3l As we will describe, the essential details in 
each step can be found in the earlier work [5l [211 125] . 

We start with the deformed conifold geometry given by the equation 

Z1Z2 = w = Z'iZi + n, Zi,weC, (2.7) 



where /z is the complex structure parameter that we take to be real positive. The geom- 
etry has the structure of T^ x M fibration over M^. Let us denote the basis cycles of T^ 
by a and /3. In the base M^, a degenerates along one line and (5 degenerates on another. 
The minimal S*^ is obtained by fibering this T^ along a line interval that connects the 
two loci. 

We wrap N branes on M = 5''^ thus engineering the U{N) Chern-Simons theory. In 
addition we place a stack of P braneqj wrapping a non-compact three cycle L of topology 
M^ X S^. The cycle L intersects M along a circle that is identified with a. These branes 
were introduced in [23] where the partition function obtained after integrating out the 
bifundamental M-L strings was shown to be a generating function for Wilson loop vevs. 
The generating function is a summation over representations of U{N) and the S^ common 
to M and L is the defining curve of the Wilson loop. Since L is non-compact one should 
enforce a boundary condition at infinity for the gauge field on the stack of branes which 
wrap L. In |23] this was implicitly done by fixing the background holonomy of the gauge 
field along a. 

A different boundary condition isolates a single Wilson loop in the representation R 
[21]. So this brane construction is equivalent to the Wilson loop insertion. See Fig- 
ure [3](a). This boundary condition is equivalent to the gauge field having a nontrivial 
holonomy matrix along the j3 cycle which encodes the data of R. 
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Figure 3: (a) The web diagram for the deformed conifold. a and /5 degenerate along the 
horizontal and vertical lines respectively. The dashed line represents S"^ that N D-branes 
wrap. The other dashed line ending on the vertical solid line represents a non-compact 
cycle L = R^ X S"^ that P non-compact D-branes wrap, (b) P non-compact branes are 
distributed along the horizontal line where a degenerates. 

First, the above brane configuration is equivalent to another system that has a new set 
of non-compact D-branes, distributed along the locus where a degenerates [25]. The new 
system has only A^ — P{= nm+i) D-branes wrapping the S^. As we review in Appendix 
[Bj a stack of rij non-compact branes sits at distance a/ = gs{Li — Lm+i) away from the 
5*^ for / = 1, . . . , m. See Figure [3](b). 

Second, by considering the new ambient geometry of Figure [3](c) with more complex 



^ Recall that P is the number of rows in R. 



structure moduli given by 

m 

ziZ2 = w, 2:32:4 = (/i-w)JJ(l -w//i/), (2.8) 

7=1 

the non-compact branes can be compactified without changing the physics. This is 
a legitimate maneuver since it reduces to the deformed conifold (12.71) by making the 
complex structure moduli ni infinite and A-model depends only on Kahler moduli. The 
result is the D-brane system from which we can derive the matrix model (12.51) . 

We now have a daisy chain of Chern-Simons theories all of them on an S^ and there are 
then annulus instantons which connect them [18] . The representation R of the Wilson 
loop determines all the necessary data, in particular the /-th Chern-Simons theory has 
gauge group U{nj), J = 1, . . . , m + 1. We get annulus instantons by integrating out the 
massive bifundamental open strings [23] . Since the mass of the string between the J-th 
and the J-th spheres is aj — aj, the interactions generated from such annulus instantons 
are summarized as 

/m+l 
Y[ [DAi]e'^^^^^'^ Y[ det (e'^^^'-^'^u} ® f/J^ - e^^^'-^^^^^f//^ ® U]) . (2.9) 
7=1 KJ 

where Scs is the Chern-Simons action, and Uj = P exp ^ Aj is the holonomy along the 
unknot in the /-th 5*^. Given this field theory description, we can now reduce it to a 
matrix model 1101: 



{Wnjr^ I IJc/HM^^e 2s» ^ -• Ijdet 2smh^ ^-^ ^ .(2.10) 

■^ 1=1 KJ ^ ^ ' 

By redefining u^^"^ — > m*^^) — QgLi = u^^^ — aj + (/-independent), we finally obtain (12. 5p . 
It is a nontrivial consistency check that the physical derivation here gives the values of 
holonomy a/ that we need to agree with the algebraic derivation in the next subsection. 
Third and finally, we can go one step further in the target space analysis though we 
have completed our task in this subsection. When each S^ in Figure M^c) undergoes a 
conifold transition the resulting closed string geometry is the toric Calabi-Yau manifold 
whose web diagram is shown in Figure [2](a). This is the Calabi-Yau manifold which is 
referred to as the bubbling geometry |5j. We will see in subsection 12.41 that the eigenvalue 
dynamics in the matrix model demonstrates the geometric transition. 



2.3 Algebraic derivation of the matrix model 

We now provide an algebraic derivation of (I2.5p . Our starting point is (12 ■4p . Using a 
standard formula for the character of U{N) this can be written as 

where Rj is as usual the number of boxes in the j-th row of R. Now we expand this ratio 
of determinants into something more compatible with the matrix model: 

, ^"^^ ,;v . ^ = 'rsgn(a)rre(^+^-^)«^«/TT(e"'-e"0 
det(e(^-j)'^0 A^ ^ "^ ' Li. ' Li.^ > 

= V" ]Te(^+^'-^)"-w/]T(e"-w -e"-W)). (2.12) 

Since Ui are dummy variables the summation over the permutation group Sn produces 
A^! identical terms, so we can write 

i i<j ^ ^ i i<j 

At this point the Wilson loop insertion has been recast into a linear term in the ex- 
ponential and a certain denominator term. There will be partial cancellation of this 
denominator term against the measure and also against the linear term. We relabel the 
variables as 

(v. V.A - (V^^^ 7/(^) 7/^^^ 7/(2) 7,(™+l) 7/('"+l)l 

{Ui, . . . , UN) — {Ui ,... ,Un^,Ui , . . . ,Un^,. . . ,u^ , . . . , M„^_^i ) 
where we recall that the Young tableau R has m blocks of rows. Then 

„ m+l ni / (I) _ (/) \ 2 / {!) _ (J'' ^ 2 

w = /nn*.'"nn 2-""^^^^ nn 2-""^^^V^ 



/=1 J=l / i<j \ / I<J i,j 

X e-2^E/..(«l")' TTTTg(JV+i^/-(JV„_,+2+i))«r' 

/ i 
/ i<j I<J i,j 



The integers Ki and Nj are defined in Appendix |Xl This can be further simphfied, using 
the trivial fact that integration variables are dummy variables, to 

{Wr) 
m+l / . / „(/) „(/)^ 2 



nh^n*."'n(2-" 

7=1 \ i i<j 



^'' ^' > e"^^'("^")'e(^+^^-^— ^+2^"^)^'"'^'^ 



'I&Snj i<j / I<J i,j 



(I) (J)\ 2 



^'^ A'^ 



(Ti€s„j Kj / Kj 1.1 e * — e J 

m+1 / , / (7) (7) \ 2 



nu7n*."'n(2-"i>- 

7=1 \ i i<j 



g 2gs ^«^ * ' e 



X 



2sinh^ 2 






1111 . „(-')^„w 

KJ i,j 



i — e J 



n Urin^^i'^n psinh^^^^-^ e-^^^^^'^^^^^--^ 

7=1 V^^' i i<j \ J 

xYlUhsinh ' ' ], (2.13) 

where Lj are defined in (12.61) and we have use the relations 

/ (.n _ (^) \ ^ 

n n 2 sinh "- ^^^- = e^^("-^)s«-''^ n n(^"''" - ^"^'v, 

7<J i,j \ / KJ i,j 

KJ i,j KJ i,j 

At this point we essentially have an m-matrix model with interactions between the ma- 
trices given by the last line in fl2.13p . 

2.4 Spectral curve as the bubbling geometry dual to a Wilson 
loop in S^ 

Matrix models have an associated geometry called the spectral curve. One can think 
of {Wr) as a single Gaussian matrix model with somewhat complicated insertion, or 
alternatively as we have demonstrated, as an ?7i-matrix model with certain simpler inter- 
actions. Taking the latter point of view, we now derive the spectral curve and explain 
its string theory interpretation. 
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M) 



The equations of motion for u] are 

(^) (I) 1 (^) (J) 

= ~uf>+gsLj + g,22^o^^ ' + :^9s J] coth ' ' ■ (2.15) 



2^ 



To solve them we define the resolventcl 

v^^\z)=g,Y.^ ' <z) = Y^v^^\z). (2.16) 



e i 



i=i e-^^ - e j=i 

We now assume that the eigenvalues distribute themselves into m distinct cuts along the 
real axis, then write (12.151) an equation on the J-th cut: 

Y.^j + Y.'^A. (2.17) 
j=i j=i / 

where v^ (z) and v_ (z) are the values of v^^\z) just above and below the cut, respec- 
tively. It will be convenient to rewrite this as 

(m m+1 \ 

Y^kj + Y,nj\. (2.18) 

j=i J=i J 

To derive the spectral curve, we generalize the complex analysis technique used in [26] 
to solve the Chern-Simons matrix model for S'^/Zp. The crucial step in solving this model 
is to find a set of functions of the resolvents v^^^ which are regular on the whole Z-plane 
where Z = e^ is of course C* valued. Then the asymptotics of v^^^ will allow us to fix these 
functions exactly and finally extract the equation for the spectral curve. The technical 
reason that we will be able to solve this model exactly is that the interaction terms in 
the equation of motion can be written polynomially in terms of the resolvents. This is 
not the case for the related A/" = 4 Yang Mills matrix models described in appendix [D] 
and also in [TT] . 

We first define some new quantities 

Xo{Z) = Ze\ 

Xi{Z) = Aie-''^'\ / = !,..., m + 1, (2.19) 

where Z = e^ and Aj = exp Qs {J2^=jkj + Yl,"j=i '^j)- Equation (12.181) implies that Xq 
and Xi are exchanged as one goes through the /-th cut, leaving any symmetric polynomial 



^ The resolvents lu^^' — gs X^i^i coth "' — in another natural definition are simply related to the 
i;(^) as w(^) =gsni -2v'^''>. 
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of {Xq, Xi, . . . , Xm+i) invariant under the process. The symmetric polynomial is regular 
on all of the cuts, and the only singularities are at Z = oo. Let us now recall the definition 
of the j'-th elementary symmetric polynomials Ej-. 

Ej{xi,...,Xn)= ^ Xi^...Xi^. (2.20) 



ji<--<« 



Together with the definition fl2.16p of the resolvents, the asymptotics as z —>■ ±00 deter- 
mine the Ej{Xq, . . . , Xjn+i) exactly in terms of Young tableau data: 



^5 



Eo{Xo, . . . , Xm+l) = 1, 

Ej{Xo, ..., Xm+i) = aj,o + aj,iZ for j = 1, . . . , m + 1, (2.21) 

Em+2{XQ, . . . , X^+i) = Ai . . . Ara+xZ. 

The coefficients are given by 

%,o = X] Bj^...Bj^ for j = l,...,m + l, 

l<Jl<--<Jj<m+l 

a,-i = Y^ Aj,...Aj^_, forj = 2,...,m + l, ai,i = 1, (2.22) 

l<,/i<-<Jj-l<m+l 

where we introduced Bi = exp Qg (X] j=/ ^J + S™=/+i '^j) ■ 

In fact the Ej appear as the coefficients of Y^ in the expansion of the function 

in+l 

f{Y,Z) ^ l[{Y-Xj{Z)) 

m+2 

= 5^(-)^r'"+2-^E,(Xo,...,x„+i) 

j=0 

m+1 

^ y-+2 + ^(_l)iy™+2-i^.^ + ^^.^^) + (_l)m2^^_4^^^^^ (2.23) 

and this vanishes upon substituting Xj for Y. So we arrive at an equation for the spectral 
curve of the matrix model (12. 5p : 

f{Y,Z) = 0, (2.24) 

where (Y, Z) are C* valued variables. 

Since f{Y, Z) is of degree m + 2 in F, the spectral curve is obtained by gluing m + 2 
cylindrical sheets. In particular (12.241) is satisfied by the total resolvent v{z) through 
substitution Y = Xq = e^'^'", and the sheet on which v{z) is naturally defined has m + 1 
cuts. By going through the J-th cut (J = l,...,m + 1), one moves to the J-th sheet as 
v{z) changes to —z — v^^'' + const. See Figure ll](a). 
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(b) 

Figure 4: (a) The spectral curve is constructed by gluing one sheet to ??i + 1 other 
sheets through m + 1 cuts. Each sheet is a cylinder parametrized by z with identification 
z ^ z -\- 2ni. Compare with Figure [2](a). (b) The vertices plot the monomials V^Z^ 
appearing in the equation (12.241) for the spectral curve. By connecting the vertices by 
suitable edges, one obtains a graph that is dual to the toric web for the bubbling geometry 
shown in Figure [2](a). 



This Riemann surface is related to a Calabi-Yau threefold in a way which is by now 
well known, namely the threefold is given by 



wx = f{Y,Z] 



(2.25) 



where w,x are C valued. It is a feature of the mirror symmetry work of Hori-Vafa 
that we can write down the toric fan directly from the Riemann surface data above. The 
recipe is to insert a vertex (a, b) on the integral 2-dimensional lattice for each monomial 
Y'^Z'' appearing in fl2.24p . By connecting the vertices with suitable edgeqj one obtains 
a graph, and the three-dimensional cone over this graph is the toric fan of the bubbling 
Calabi-Yau. The two dimensional graph is the dual graph of the toric web diagram, so 
from Figure m^b) we see agreement with the previous work [5]. 

For concreteness we now work out the simplest case when i? is a rectangle. In this 
case the nontrivial data is 









(2.26) 



^ In the limit of large Qsni and Qskj (the large volume limit in the S^ case, but not in the S'^ j'L^ 
case), the difference between the GLSM algebraic coordinates [28] and our moduli is suppressed. The 
mirror curve of the GLSM in this limit agrees with our spectral curve including the coefhcients, with 
the choice of internal edges in Figure IHb) 
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and 



ai,o = l + e*^+*% ai,i = l, (2.27) 

02,0 = e*^+*^ a2,i = e*i+*^+*^ + e*^ 
a3,o = 0, a3,i = e*^+*^+2t3^ 

with ti = Qsni, ^2 = S's^i, ^3 = 5's''^2 5 and so the spectral curve is exphcitly given by 

r^ - (1 + e*2+*3)F2 _ y2^ ^ gt2+t3y ^ (^gii+t2+t3 ^ e*^)rZ - e*i+*2+2t3^ ^ q^ (^2.28) 

2.5 Eigenvalue distribution 

The exact eigenvalue distribution can be obtained by solving fl2.24p for v{z) via Y = 
exp(z + v) and by computing the eigenvalue density p oc v+{z) — V-{z) along the cuts. 
Here we apply force balance to derive the approximate distribution when 

g^rij > 1, gskj > 1 for all /. (2.29) 

Force balance is easier to understand intuitively. 

We make the assumption, to be justified a posteriori, that 

uf ^ - nj^^ > 1 for all J, J, i,j such that I < J. (2.30) 

Because the last term in (12.151) becomes constant we have 

n I m+1 



'^^" =9sY,- MT-m + 9s[Y,kj -Y.'^J + Y. ""A- (2.31) 



1-e^- -"» 



2 

j^i I — C 3 ' \J=1 J=l J=I + 1 

We expect that when ggUj is large, the eigenvalues of u^^^ spread over a large region 
allowing us to approximate the function 1/(1 — e^) in (I2.3ip by a step function. If we 

/ m+l 

M) 

— ^ysi' -r ys i 

vJ=l J=l J=/+l 



order the eigenvalues so that u\ < u\ for any i < j, it follows that 



Ui 



2gsi + gs\^kj -^nj+ ^ n,j\ , i = l,...,n/. (2.32) 



Along the /-th cut that has width 2gsni, the eigenvalues of m*-^^ are distributed uniformly. 
The J-th and / + 1 cuts are distance ggkj apart from each other We can thus justify 
the approximations above when ggfij and gskj are all large. See Figure [51 As discussed 
above, this sheet is connected to other m + l sheets through the m + l cuts as shown in 
Figure 111(a). 



^ Since u is the holonomy along a, its eigenvalue distribution is different from the distribution (JB.63P 
of holonomy ^„ A. In particular the eigenvalues are quantized in unit of 2gs . It should be possible to 
physically explain (|2.32p using the fact that the matrix model captures the Wilson loop in a non-canonical 
framing |10j. 
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Figure 5: The eigenvalues are distributed along m + 1 cuts on the cylinder parametrized 

by z. 



3 Bubbling Calabi-Yau for lens space from a matrix 
model 

3.1 Matrix model for a Wilson loop in lens space 

A simple generalization of the topological A-model on T'S'^ is the orbifold Xy = T*{S^/Zp) 
[TU] . The particular orbifold action is such that S^/Zp is the lens space L{p, 1). This 
space is defined by the equation 

|zi|2 + |z2r = l (3.33) 

for complex variables Zi and Z2, together with identification 

(zi,^2)~(e''^'/^zi,e-2-/Pz2). (3.34) 

We study the Wilson loop 

Wr = TiRPe^^ (3.35) 

along a circle that is the generator of the fundamental group. We assume that the circle 
is the unknot. 

The U{N) Chern-Simons theory on L{p, 1) has many vacua. Since the equation of 
motion is solved by a fiat connection, the vacua are in one-to-one correspondence with 
the A^-dimensional representations of 7ri{S^/Zp) = Zp. The group Zp is abelian, so any 
such representation is a sum of one- dimensional ones. A one- dimensional representation 
is specified by an integer a = 1, . . . ,p. Thus a vacuum is specified by a partition of A^: 

N = Ni + N2 + ... + Np. (3.36) 

Here Na is the number of times the a-th irrep appears. The contribution of this vacuum 
to the partition function is given by 

Zp = lfldu,ll (2 sinh ^^) ' exp f-^ E «' + ^ E ^^"^) ' (^.37) 
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This matrix model was formulated in [9] , and was studied for example in 

According to the prescription in [10] (see also [31]), the contribution from this vacuum 
to the Wilson loop vev is given by 



i=l i<j 

where n is a vector of integers 



{WR)p=ff[du,l[(2smh^^^^j exp I -^Y^u] + — J^n^Mi j Trfidiag(e"0, 

(3.38) 



n = (C^,C^, • • • ,^7"^). (3.39) 

The spectral curve for this matrix model can be derived [26] and it agrees with the string 
theory prediction [29] . 

For a large representation R with large values of nj and kj, we expect a large backreac- 
tion of fields to the Wilson loop insertion. We propose that the gauge field path-integral 
has now more saddle points. Each saddle point specified by (A''^) before insertion splits 
into many each of which is specified by non-negative integers {Nia) satisfying the con- 
straints 

p m+l 

Y, Nia = ni, Y ^1- = ^- (3-40) 

a=l 1=1 

We will argue that the contribution to the Wilson loop vev from the saddle point 
specified by the (iV/a) is given by the multi-matrix model 



{Wr)\: 

. m+l p 



(Nla) 



I n n^H«(^")exp (-^T,[u^^-)f + (lj + —a] Tiu^'A (3.41) 

•^ 1=1 a=l \ 9s \ Us / / 

JJ det (2 sinh """^l -!»'■'"') n de' (^ -•> """ "' ^ , ' ^ "''") ^ • 



I,a<b ^ ' I<.J,a,b 

Wilson loops in the lens space matrix model have also been considered in the interesting 
recent work [31] and it would be of interest to apply their methods to the spectral curve 
in this paper. 

3.2 Physical derivation of the matrix model 

We now derive the matrix model from a D-brane configuration that realizes the Wilson 
loop in a lens space. 
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Let us recall that Xp is a Zp orbifold of the deformed conifold given by (12. 7p . The 
orbifold action is generated by 



(e-2"*/Pzi, e2"^/P;Z2, e^'^'^^zs, e^'^'^'l^'z^) 



[zi, Z2, Z3, Z4) ^ [e "■""'^1,6""'' ^2:2,6'-""^ 2:3,6 '■"•■"'Z4), (3.42) 

and the Zp action on the S^ given by ^2 = ^1, -^4 = —-^s (so |2;3p + |2;ip = yu) defines the 
lens space L{p, 1) = S'^/Zp. Since the Zp only acts on the phases, Xp is still a fibration 
of T^ X M over M^. Let us redefine a to be the 1-cycle corresponding to the generator 
of the fundamental group, and /? the 1-cycle given by the 27r phase rotation of z^. We 
use the axes of the two cylinders (given by Z1Z2 = const., ^3^4 = const.) and the Re{w) 
direction as the base M^. The cycle /3 degenerates at w = fi and so does /3' = —pa + [3 
at w = 0. The cycle a never degenerates. 




i^m+l. 



nm:r ^2/V y/V rii 




(a) (b) (c) 

Figure 6: (a) The cycle /3 degenerates along the vertical line while —pa + (3 degenerates 
along the other line. If a linear combination qa + r(3 degenerates, it does so along a line 
in the (g, r) direction, (b) P non-compact branes are distributed along the line where 
—pa + [3 degenerates, (c) There are m + 1 copies of S'^ /Zp. 



We engineer U{N) Chern-Simons theory by wrapping A^ D-branes on the S^/Zp. To 
insert a Wilson loop along the knot a, we consider P D-branes that wrap the non-compact 
cycle L = M^ X S*^ in which (3 is contractible. See Figure [n]( a) . The boundary condition 
{R\ on the P branes picks out the Wilson loop insertion in representation R, as explained 
in [5]. The boundary condition induces holonomy 



/3=pa+/3' 



A = diag lgjRi-i + -{P + N + 1) 



(3.43) 



i=l 



along the contractible cycle (3 = pa + (3' . By fibering the T^ over a semi-infinite line 
ending on the locus where (3' degenerates, we obtain a 3-manifold in which (3 is non- 
contractible. We can consider a configuration of D-branes wrapping this 3-manifold. Is 
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the configuration equivalent to the one we started with, as in the S^ case? We assume 
it is, and we will see evidence below. The basic nontrivial cycle in the new 3-manifold is 
a, and the holonomy along it is given by 

A = idiag (gs(R^-^ + ^{P + N+ 1) j j (3.44) 

because /5' is contractible. See Figure [6t^b). 

As in the S^ case, it is natural to split the P non-compact branes into m stacks with the 
J-th stack containing nj branes. We can now replace Xp by the Zp orbifold of the large N 
dual geometry given by the equations (12.81) . This is possible because (12.81) are invariant 
under the orbifold action. The non-compact branes are now replaced by compact ones 
wrapping copies of lens space S^/Zp. Thus we reach the desired system of D-branes, 
whose world-volume theory is m + 1 copies of Chern-Simons theory on lens space S^/Zp, 
interacting via Ooguri-Vafa operators. The system is shown in Figure [6](c). 

To write down the matrix model, we need to choose the vacuum of the theory. We 
have a U{ni) Chern-Simons theory on the /-th lens space. As reviewed in the previous 
subsection, the theory has many vacua corresponding to the choice of a fiat connection. 
Let us choose the vacuum specified by the partition rii = ^^ Nja. Then according to the 
prescriptions in [lOj, the contribution to the Wilson loop vev from this vacuum is given 
by 

"^ I,a,i I,a<b ^ 

TT det (2 sinh (^^^"^ + (^iIp) ®^-^® (^^^''^ + Qj/p) 

^ 1 n V, \ 



X 

I<J,aj 



X exp f-^Tr(u(^"¥ + —aTiu^'A . (3.45) 

V 2^, gs J 

By redefining the variables as u^^""^ -^ u^^""^ — gsLj/p = u^^""^ — aj/p + (/-independent), 
we obtain (13.411) . It is remarkable that we get the holonomy aj/p, including the factor of 
1/p, which is necessary to be consistent with the algebraic derivation. The success gives 
us confidence in the assumption we made above. 

This brane construction makes it clear what the dual bubbling geometry should be. It 
should be the toric Calabi-Yau shown in Figure [2](b), where all the copies of lens space 
have undergone geometric transition. This proposal will be confirmed in subsection 13.41 
by deriving the spectral curve of the matrix model, and showing that it is the mirror of 
the toric Calabi-Yau. 
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3.3 Algebraic derivation of the matrix model 

The vector of integers n in (13.371) breaks the U{N) invariance down to the product 
subgroup XaU{Na) and subsequently the Sn symmetry to S' = XaSNa- Nonetheless the 
Wilson loop is in the representation R of U{N) and as such we cannot immediately apply 
all the steps we used to solve the S^ case in section 12.31 The workaround is to consider a 
generating function of matrix integrals, one term of which will correspond to the Wilson 
loop vev {Wji)p ''"" . This generating function will have S^ symmetry and thus we need 
only the technology used in section 12.31 to solve this case as well. 
So we will consider the generating function with variables zi^ . . . ,Zp 

WR^p{Za) 
■^ ' i=l i<j ^ ^ \i=l a=l J 

The coefficient of Yla ^a"" i^ (13.461) is (yVii)p. Since all the Ui are dummy variables on the 
same footing, we can straightforwardly repeat the analysis of section 12.31 to arrive at 

„ m+l ( . ni / (/) _ (/) \ 2 

X n n (e ""'•" -) n n (^ >""'> ^^^^^) . (^ad 

I i \ a J KJ i,j \ J 

where the eigenvalues (ui) have been divided into m groups (ul ) (?7i is again the number 
of groups of rows in R). To understand the coefficient of H-^a^" i^ i^ best to divide up 
the eigenvalues (mj) into {ul ) and also {u] ) such that 

p m+l 

M)\ _ I l^„,(^»)^ ^„,(«)^ _ I \/„M^)\ 



U(-!"^ (-S"^)=U 



[Ui ) = \_\[ul 'j, [Ul ') = U[ul ,, 

a=l 1=1 

m+l p 

in.) = L|K^'^) = U(^S"^)- (3.48) 

1=1 a=l 

So clearly we have the constraints (13.401) and for each choice of non-negative integers 
{Nia) which satisfies these constraints, we have the following contribution to {Wr)p. 

(M-'r)'"'-' 



n 



a,i a,i<j \ / a<b,i,j \ / 



X exp -— 2_^{u\ ') + 2_^{Lj + a)ul 

9s ■ . 9s 
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Yl 2 sinh 

I<J,a,b,i,j 



(la) (Jb) - 

U\ - Uj 



This is the matrix model fl3.4ip in the eigenvalue basis. We now solve the matrix model 
and derive its spectral curve. 

3.4 Spectral curve as bubbling geometry for a Wilson loop in 
lens space 

We now derive the spectral curve associated to f l3.49p that captures the contribution of 
the particular vacuum specified by the integers (Nia)- Since the gauge theory sums up 
such contributions, the Wilson loop is actually dual to a sum over geometries. 
The equation of motion for u- that follows from fl3.49p is 

(la) (la) (la) (lb) 

(I a) v^ , K '^j sr^ . K '^i 

= -pul ' + gsLi + 2-nia + Qs 2_^ coth — '- \- 9s 2_^ coth — '- — 

j^i bjta,i 

+ -(7,J] coth^ ^-^, (3.50) 

J^I,b,i 

and so we first define several resolvents 

Nia „(-?") P m+l 

j=l 6 ' 6^ a=l 1=1 



In terms of these we can write ( I3.50p as an equation on the (Ja)-cut: 

(m m+l \ 

^ fcj + ^ nj J + 27rm. (3.52) 

j=i J=i / 

Following the same procedure as in section 12.41 we define some new variableqj 



Xo = ZPe" 



Xj = Aje-""", J = l...,m + 1, (3.53) 



where Z = e^, Ai = exp gg {Yiy=i ^J ^ Y17=i "^j)- Then the spectral curve is agai 



am 



given by (recalling once more that (Y, Z) are C* valued variables) 

/(r,Z)=0 (3.54) 

where 



m+2 



/(F,Xo,...,X„+i) = \{{Y-X,) 

m+2 

= ^(_).y™+2-,i5;^.(Xo,...,X„+i)- (3.55) 



j=Q 



^Despite identical nomenclature these variables are of course unrelated to those in section \2M 

19 



The difference with section [23] hes in the asyniptotics of the elementary symmetric poly- 
nomials, from which we can determine their structure: 



i=0 



Eo{Xo, . . . ,Xm+l, 
Ej{Xq, . . . ,Xm+l, 
Em+2{Xo, . . . , Xm+l) = Ai . . . Am+l 

Some coefficients are easily determined: 



aj^iZ' for J = 1, 
ZP. 



,m + 1, 



(3.56) 



*j>p 



J2 Bj,...Bj^ forj = l 
i<Ji<-<j,_i 



Aj. 



j-i' 



. . ,m + 1, 
for j = 2,...,m + 1, ai,p 



where Bj = expgs {Yl^=i ^J ^ jyj=i+i^j) ■ The remaining aj^i are complex structure 
parameters that are determined by demanding that f v{z)dz = —2TTigsNja for the integral 
around the (Ja)-cut. We can again write down the toric fan of the bubbling Calabi-Yau 
geometry directly from the spectral curve by plotting the monomials V^Z^. See Figure 
IHa). 





(b) 



Figure 7: (a) The vertices plot the monomials Y"'Z'' in ( I3.54p . By connecting the vertices 
by suitable edges, one obtains a graph that is dual to the toric web in Figure [2](b). (b) 
The eigenvalue distribution on the cylinder. Here we chose p = 3 for illustration. 



We see that this toric threefold is a daisy chain of lens spaces, and the role of the 
complex structure deformations of the spectral curve is to desingularize each lens space. 
An interesting new feature of this geometry is the presence of nontrivial four-cycles. 



3.5 Eigenvalue distribution 

As in the S^ case, when all gsUj and gskj are large, the interactions between u'^^^ and 



u 



(J) 



can be neglected. The eigenvalue distribution for each / is then that of a single lens 
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space obtained in [26]. According to [26], the eigenvalues of m^^"^ are distributed along a 
cut at Ini(z) = 2TTa/p parallel to the real axis. See Figure [Tl^b). This sheet is connected 
to TTi + 1 other sheets, each through p cuts. The resulting topology is that obtained by 
fattening the toric web in Figure [2](b). 
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Appendix 

A Summary of Young tableau data 

The Young tableau R has n/ rows of length Kj such that Ki > K2 > . . . > K^ > 

Km+i = 0. It also has kj columns of length Nm-i+i such that Ni > N2 > . . . > N^- We 
also define n^+i = N — J27=i ^h^o — ^5 ^^^ Km+i = 0. The integers ni,ki,Ni, and 
Kj satisfy the relations 

m-I+l 

Ni= Y^ nj for J = 0,1,..., m, (A.57) 

j=i 

and 

m 

Ki = Y,kj for / = l,2,...,m, K^+i = 0. (A.58) 

j=i 

See also Figure [H We also denote by P the number of rows in R, so P = Ni. 
Other useful sets of quantities are 

m -. /—I ^ m+1 

j=i J=i J=i+i 

ai = QsIkj- lni + ... + nj^i + -nj\ +-{P + N) 

= gs{Li-Lm+i), (A.60) 
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and 

n m+1 



Aj = exp(7J^fcj+^nj , (A.61) 

\j=i j=i J 

/ m m+1 \ 

Bi = exp(7J^/cj+ 5Z rij . (A.62) 



yj=l J=I+1 



B Area of the annulus diagrams 

Here we explain the identification aj = gs{Li — Lm+i) in subsection 12.21 

The P non-compact D-branes with the boundary condition {R\ has background holon- 
omy [25j (gauge equivalent to the position relative to the N — P compact branes on 

j>A = diag {gs{R^ -i+'^{P + N + l))\ (B.63) 

along the /3 cycle. When we split the P non-compact branes into m stacks, the average 
value of the holonomy in the J-th stack is 

ai = gAKi-Ui + ...+ni^i + -nA+-{P + N)\, I = l,...,m. (B.64) 

Since this is the distance from the S^, it is natural to define a^+i = 0. The parameters 
ai (/ = 1, . . . , m -|- 1) are then the positions of m + 1 copies of S^ in the new geometry 
given by (12.81) . aj — Om+i is the area of the annulus between the 5*^ and the /-th stack 



of non-compact branes. See Figure [3]^b). Note that (IB.64P can be written as aj = 
gs{Li - Lm+i) = 5's-^/+ (/-independent). 

C Alternative matrix models for a Wilson loop in S^ 

Here we discuss two alternative matrix models whose partition functions are the Wilson 
loop vev for Chern-Simons on S^. The first is 

(Wr) 

duu dU^^UU^^^ . . . dU^"'^e-^^'^'''\det U^^^f-{det f/(2))fc— i . . . (det U^""^)^' 



^ det(l - e« ® t/(i)-i) det(l - f/(i) ® t/(2)-i) ' ' ' det(l - f/(™-i) ® f/M-i) ' ^^'^^^ 
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Here t/'-^-* is an Nj x Nj unitary matrix. The second is 

(Wr) 

X det(l + e" « C;<..-.)^-^^-^A___ . . . ___l___^c.66) 

for which U^^^ is a Kj x i^/ unitary matrix. 

These models are obtained from (12. 4p by the same algebraic manipulations that led to 
similar multi-matrix models for A/" = 4 Yang-Mills in [11] . 

C.l Physical derivation 



Here we give a physical derivation of the matrix model (1C.65I) from a D-brane configu- 
ration. 

We begin with the configuration of A^ compact and P = Ni non-compact D-branes 
(Figure [3](a)) that we discussed in subsection (12.21) . On the non-compact branes we 
impose the boundary condition {R\ to picks out the Wilson loop Wr from the annulus 
diagrams between the branes. 

We now consider a new configuration that realizes the Wilson loop insertion. We mod- 
ify the geometry and introduce another locus on which (3 degenerates. By fibering the T^ 
over a line interval that connects the two loci where (3 degenerates, we get a cycle of topol- 
ogy 5*^ X 5*^. We wrap A^^i D-branes around this cycle while placing external fundamental 
strings in an appropriate configuration. This configuration of the fundamental strings 
is that they insert the Wilson loop in the one-dimensional representation A^^™ j25jJ1 
Additionally we place N2 non-compact D-branes that end on the second locus where f3 
shrinks. We choose the boundary condition to be (Q*-^-*], where the Young tableau Q*-^^ is 
obtained from R by removing the first k^ columns (Figure [9]). The external strings and 
annulus diagrams from the non-compact branes insert to the 5*^ x S"^ branes the Wilson 
loop 

Tr^».™e^^TrQ(2)Pe^^ = Tr^Pe^^. (C.67) 

Since S^ x S*^ is obtained by gluing two copies of solid torus by identifying their bound- 
aries, the path- integral there reduces to the inner product. Thus from the annulus di- 
agrams between the S^ and S^ x S"^, the path- integral picks out the combination that 
inserts the Wilson loop Wr into S^. See Figure [8](a). 



^ Ani is the rank A'^i totally anti-symmetric representation of U{Ni) and is one-dimensional. 
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N, 



Ni 



-W^) 



No 



N, 



Ni 



INo 



(a) 



(b) 



^Nr, 



Figure 8: (a) The P = Ni non-compact D-branes in Figure [3] (a) are compactified by 
modifying the Calabi-Yau geometry without changing the topological string amplitudes. 
The state \R) specifying the boundary condition is implemented by placing external 
string world-sheets that insert the Wilson loop Tr^j exp ^ A. (b) The geometry and 
the configuration of D-branes and non-compact string world-sheets that give rise to the 
multi-matrix model (1C.65I) . Each horizontal dashed line represents D-branes wrapping 
a Lagrangian sub manifold of topology S^ x S"^. The cylinder ending on the J-th dashed 
horizontal line represents fundamental strings in a configuration that inserts a Wilson 



loop in the representation A 



N, 



-l+i 



for / 



,m. 



We can repeat this process (Figure [9]) and show that the following configuration is 
equivalent to the Wilson loop insertion. The total geometry is given by the same equation 
(12. 8p as in subsection 12.21 with one locus where a shrinks, and m + 1 parallel loci where 
(3 shrinks. A^ D-branes wrap the original S^. We also wrap Nj D-branes on the S^ x S"^ 
between the J-th and (/ -|- l)-th loci where j3 shrinks. Finally we place fundamental 
strings, along the J-th locus, that insert the Wilson loop in the representation Aj^^""'' 
into the J-th S^ x S'^. See Figure [8](b). 




iVi 



A^2 




■N, 



m—l 1 



NmO 



W 



R = Q^^^ Q(2) . . . Q(m-1) Q{m) 

Figure 9: A shrinking sequence of Young tableaux R = Q*^^^ D Q^^-* 'D . . . ^ Q^"^\ 



Using the prescriptions in [lOj, we obtain the matrix model (JC.651) from this D-brane 
configuration. There is no Gaussian factor for the Chern-Simons on S^ x S"^ since the 
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path- integral is simply the inner product. The external fundamental strings insert the 
determinant factors. 



It is also easy to extend the derivation to ( 1C.66I) . this time using anti-branes instead 
of D-branes. This explains the appearance of one determinant, rather than the inverse 

of it, in flaeej) . 



C.2 Solving ( ICTMI ) 



Now that we know the physical origin of the matrix model flC.651) . let us here solve it in 
the large N limit. In terms of the eigenvalues, the matrix model can be written aqj 



m Ni 



/lY III, -^ * J -1 -i • '"' -''J 

+ XI ^°S f sinh ' ^ j + X] 5Z ^°S ( sinh 



u. 



(7) (7)x 2 



(I) 



2 / -^^ -^^ \ 2 

i<j ^ ^ 7=1 i<j 

N Ni m-1 Nj Ni+i 

E E i°g(i - ^"'"^'') - E E E i°g(i - ^ 

j=l j=l 7=1 i=l i=l 



(-f) (-f+i) 



(C.68) 



Proceeding as in subsection 12. 4^ by defining the resolvents 

i=\ 

Nr 









for J = 1, ... ,771, 



we express the saddle point equations as 



v±{z) + z = -v^{z) + V 



(1) 



fl's^m+1 



on the M-cuts and 



- v^'-'\z) + 7;£)(z) = -v^4\z) + v^'+'\z) + ^.(fc„,_7+i + n™_7+i) 



(C.69) 



(C.70) 



(C.71) 



on the M^'^-'-cuts, for / = 1, . . . , ?7i. Note that we have defined v^^^ = v, ti(™+^) = 0. These 
equations state that the following quantities are permuted as one goes through a cut: 



Xo = e^+^ Xj^Aje-'''""'^''^'''"'-'^'' for / = 1, . . . ,m + 1, 



(C.72) 



^ The quantities u^ and v^^' in this subsection are not to be confused with the quantities denoted 
by the same symbols in other parts of the paper. 
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where Aj is the famihar quantity defined in (lA.Gip . The asymptotic behavior of Xj as 
z -^ ±cxD is the same as that of Xj in subsection I2.4[ The rest of the analysis then goes 
exactly in the same way, leading to the spectral curve fl2.24p . In particular Xj here can 
be identified with the quantity denoted by the same symbol there. It was found there 
that Xq has m + 1 branch cuts, while Xj with I = 1, . . . ,m + 1 shares with Xq just the 
/-th cut. One can now show using (1C.72|) that v^'^~^^^\z) shares with v{z) the first / of 
these cuts, and thus the /-th cut consists of the eigenvalues of u, u^^\. . . , and 'u(™^^+^\ 
How do we interpret the different kinds of eigenvalues that lie along the same cut? We 
believe that these eigenvalues form bound states due to attractive forces, as explained in 
pT] for a matrix model that describes a Wilson loop in A/" = 4 super Yang-Mills. The 
/-th cut has rij (u-u^^^-. . .-u^"^~^~^^^) bound states. 



C.3 Solving ( ICTMI ) 

Let us also solve (1C.66I) , which in terms of eigenvalues readqj 

/N m Ki r N m Ki 

i=l 7=1 i=l L ^* i=l 7=1 i=l 



U 



in 



+ y^ log I sinh 



u~ — u. 



E E ^°S sinh 



«P-^^)~ 



i<j ^ ' 7=1 i<j 

N Ki m-1 Ki Kj+i 

+ E E i°g(i - ^"""^'') - E E E i°g(i - ^" 

j=l j=l 7=1 j=l j=l 



(I) (^+1) 



Again by defining the resolvents 

v(z) = 



N 

i=l 
Ki 



oUi 



,m 



'■'"('-) = 9.E^ 



for /= 1, ... ,771, 



i=l 



e i — e 
the saddle point equations can be written as 

{v{z) + z)^ = {-v{z) - v^^\z) + g,{N + K^))^ 

on the M-cuts, 



,(1) 



(_^(^) _ v^^\z))^ = {v^^\z) - v^^\z) - gs{ni + k,)) 



T 



(C.73) 



(C.74) 



(C.75) 



(C.76) 



^ The quantities ui and v^^ jn ^i^ig subsection are not to be confused with the quantities denoted 
by the same symbols in other parts of the paper. 
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on the ■u'^^^-cuts, and 

{v^'-'\z) - v('\z))^ = {v^'\z) - v^'+'\z) - g,{nj + kj))^ (C.77) 

on the M^'^-'-cuts for / = 2, . . . , ?7i, where we defined t>(™'+^) = 0. From these equations we 
see that the following quantities are permuted as one goes through a cut: 

X^ = e''+^ X[ = Aie"''-^*'' , X\ = Aje''^'''^-"^'^ for / = 2, . . . , m + 1, (C.78) 



where Aj are defined in (lA.6ip . The asymptotic behavior of Xj as z ^ +oo is that of Xj, 



but as z ^ — oo, X[ behaves like Xm+i, and Xj hke X/_i for J = 2, . . . , m + 1. Xq and 
Xq share the same asymptotics, hence so do Ej{Xq, . . . ,X!^_^_i) and Ej{Xo, . . . ,Xm+i)- 
One concludes that the spectral curve of this model is the one found in subsection I2.4[ 

What is the explanation of the difference between Xj and Xjl The functions Xj are 
all holomorphic on the zero-th sheet except on the m + 1 cuts along the real axis. While 

(Xq, X-^^, . . . , X^_^^) = {Xq, Xi, . . . , Xm+i) (C.79) 

for Iie{z) that is positively large enough, for negatively large Re(2;) we have 

(Xq, X^, X2, . . . , X^_^_{) = (Xo, Xm+i, Xi, . . . , Xm). (C.80) 

Thus Xj are not continuous, and we believe that the discontinuities arise due to the 
f '^'^-'-cuts (I = 1, . . . ,m) that lie in the imaginary direction as in [TT]. 

D An improved matrix model for A/^ = 4 Yang-Mills 

This appendix is targeted at readers who are interested in Wilson loops in the AdS/CFT 
context. 

It is believed [HI [32] that the correlation functions of circular loops in A/" = 4 Yang- 
Mills are captured by the Gaussian matrix model. The precise correspondence states in 
particular that 

Tr^Pexp ([{A + 9'X'ds)\ = j f dM exp (-'^TtmA Ttrc^^ . (D.81) 

The left-hand side is the normalized expectation value of the circular supersymmetric 
Wilson loop in the Yang-Mills with gauge group U{N). The right-hand side is normalized 
by using the partition function Z which is the integral without the insertion of Tr^je*^. 
dM is the standard hermitian matrix measure, and A = Qym^ is the 't Hooft coupling. 
In the absence of operator insertions, the eigenvalues are distributed according to the 
Wigner semi-circle law in the large N limit. 
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By applying the same algebraic manipulation as we did in subsection I2.3[ we conclude 
that the vev of a circular Wilson loop is given by several Gaussian matrix integrals 
correlated by interactions: 






1=1 KJ 

3+1 / T ni 

1=1 \ "^ i=l l<i<j<ni 



^ / n h^ n d"^^ n ("^^^'^ - mf ^^^.-^ ^.k^'')%^. e. -^^ 



X n nn . „».i... - i^^-^^) 

l</<J<g+l i=l j=l 



l_e"^. -™' 



Here M*^^) is an rij x n/ hermitian matrix. This is the direct analog of the second 
expression in (12.131) . We used the symbol g to denote the number of blocks in R as in 
[H [TT] , so g = m in the notation of Figure [H 

Using this multi-matrix model, it is remarkably easy to obtain the eigenvalue distribu- 
tion and reproduce the Wilson loop vevs for the representations R that are realized by 
a D3-brane p^[T3] . D5-brane [121 [H], ^.nd bubbling geometry [11]. In particular, for an 
R with large gyM^i ^^^ 9ym^Ii the gravitational dual is a smooth bubbling geometry. 
Since m\ is pulled to the right by the linear potential in (ID. 821) with coefficient i^/, 
is much larger than m) ii I < J. Then the interaction between M^^'> and M^^^'> can 



m. 



(^) ; 



be neglected. It then follows that for each M*^^) the eigenvalues are distributed around 
XKj/AN according to the semicircle law with half width \/g^M^i- 
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